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Abstract 

Physical layer multicasting with opportunistic user selection (OUS) is examined in this work for 
multicell multi-antenna wireless systems. In multicast applications, a common message is to be sent by 
the base stations to all users in a multicast group. By adopting a two-layer encoding scheme, a rate- 
adaptive channel code is applied in each fading block to enable successful decoding by a chosen subset 
of users (which varies over different blocks) and an application layer erasure code is employed across 
multiple blocks to ensure that every user is able to recover the message after decoding successfully in a 
sufficient number of blocks. The transmit signal and code-rate in each block determine opportunistically 
the subset of users that are able to successfully decode and can be chosen to maximize the long-term 
multicast efficiency. The employment of OUS not only helps avoid rate-limitations caused by the user 
with the worst channel, but also helps coordinate interference among different cells and multicast groups. 
In this work, efficient algorithms are proposed for the design of the transmit covariance matrices, the 
physical layer code-rates, and the target user subsets in each block. In the single group scenario, the 
system parameters are determined by maximizing the group-rate, which is defined as the physical layer 
code-rate times the fraction of users that can successfully decode in each block. In the multi-group 
scenario, the system parameters are determined by considering a group-rate balancing optimization 
problem (i.e., a max-min weighted group-rate problem), which is solved using a successive convex 
approximation (SCA) approach. To further reduce the feedback overhead, we also consider the case 
where only part of the users feed back their channel vectors in each block and propose a design 
based on the balancing of the expected group-rates. In addition to applying SCA, a sample average 
approximation technique is also introduced to handle the probabilistic terms that arise in this problem. 
The effectiveness of the proposed schemes is demonstrated through computer simulations. 
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I. Introduction 

A. Motivation and Background 

Multicasting has attracted much attention in recent years due to the inereasing demand for 
mass content distribution, such as software and firmware updates, file downloads, and multimedia 
streaming. In these applications, common information is to be disseminated effieiently to all users 
in a multieast group. Due to its importanee in mobile applieations, sueh service is also being 
introdueed into eurrent and next-generation eellular standards, sueh as the global system for 
mobile eommunications (GSM), the worldwide interoperability for mierowave aeeess (WiMAX), 
the long term evolution (LTE) ete., in the form of the so-ealled multimedia broadeast/multicast 
service (MBMS) [l]-[3]. Different from wireline networks, multicasting in wireless systems 
enjoys the so-ealled wireless broadeast advantage (namely, the advantage that all users within 
the transmission range ean reeeive) and, therefore, many researeh efforts have been devoted to 
the development of physical layer teehniques that ean fully exploit these advantages. 

Most works in the literature on wireless physieal layer multieasting eonsider the transmit signal 
design at the base station (BS) or BSs for effieient delivery of eommon information to all users in 
a multicast group. Due to fading, the ehannels experienced by users in the multicast group may 
vary drastically and, thus, the rate of the ehannel eode must be low enough to ensure that all users 
in the group ean suceessfully deeode. When multiple antennas are available at the transmitter, 
beamforming and precoding techniques can be further employed to improve the effeetive ehannel 
quality of the worst user [4]-[6] and, thus, the multieast rate. With no restrietions on the rank 
of the signal eovarianee matrix, the optimal preeoder (or the transmit eovarianee matrix) ean 
be found using semi-definite programming (SDP) teehniques [7]. However, to reduee decoding 
complexity at the reeeivers, many works, such as [4] and [5], considered instead the use of 
multicast beamforming (i.e., signals with rank-1 eovarianee matrices) and proposed approximate 
algorithms for the design of multicast beamformers, which is otherwise known to be NP-hard 
[4]. These works were extended to systems with multiple multieast groups in [8]-[12], where 
eoehannel interference between different groups’ messages was taken into eonsideration, and 
were also extended to multicell systems in [13], [14], where eoordinated transmissions among 
BSs were employed. However, even with the above signal designs, the effieieney of the physieal- 
layer multieast transmission is still fundamentally limited by the user with the worst ehannel 
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Fig. 1. Illustration of the opportunistic multicast scheduling (OMS) scheme. 


and the rate required for all users to deeode may eventually go to zero as the number of users 
in the group inereases. 

B. Related Works on Opportunistic Multicasting 

For delay-tolerant applications, such as file downloads and software or firmware updates, the 
limitations caused by the worst user in the network can be overcome by dividing the transmission 
into multiple blocks and by scheduling only a subset of users to receive in each block. This 
technique was referred to as opportunistic multicast scheduling (OMS) in [15]-[18]. In particular, 
the OMS scheme requires a two-layer encoding scheme, as illustrated in Fig. 1, which consists 
of a physical layer channel code (hereby referred to as the inner code) whose rate is adapted to 
the chosen user subset in each block and an application layer erasure code (referred to as the 
outer code) that is performed across multiple blocks. By selecting only a subset of users to serve 
in each block, the rate of the inner code would be less restrictive since successful decoding 
needs to be guaranteed for fewer users in each block. The overall group-rate, defined as the 
code-rate multiplied by the fraction of users in the group that successfully decode, can thus be 
effectively increased. However, since a user may not be able to successfully decode in every 
block, the channel that it experiences across multiple blocks is effectively an erasure channel, 
and an outer erasure code should be employed to ensure that all users in the multicast group 
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can eventually reeover the message after deeoding over a suffieient number of bloeks. The OMS 
seheme ean aehieve signifieant performanee gains, but eomes at the expense of delay. That is, 
these gains result from using time as an additional resouree or degree of freedom as eompared 
to eonventional sehemes that require sueeessful and instantaneous deeoding in every bloek. 

In praetiee, the outer eode ean be implemented using, e.g., LT, Raptor, and Fountain eodes 
[19], [20]. With OMS, the group-rate is known to eonverge to a non-zero eonstant as the number 
of users ineeases [15]-[17]. This is in eontrast to oases without OMS where the group-rate is 
known to diminish to zero. Most works in the literature on OMS foeused on the single-eell single¬ 
antenna soenario, as in [15]-[18]. Our work foouses instead on the multi-antenna soenario, where 
the problem is oonsiderably more diffieult due to the dependenee between the user seleetion 
and the transmit signal design. The multi-antenna soenario was also examined more reeently 
in [21]-[23], but only for single-eell soenarios. In partioular, in [21], the joint preooding and 
user seleetion was performed using a heuristio semi-orthogonal veotor seleetion algorithm. Their 
proposed seheme has low eomplexity but is suitable only for single-eell soenarios with sum 
power oonstraints. In [22], [23], the transmit signal design was restrioted to beamforming and 
the user seleetion was performed using a heuristio subset seareh algorithm. The algorithm ean 
be extended to multioell soenarios, but is subjeot to high eomputational eomplexity. 

C. Main Contributions 

The main eontribution of this work is the development of effioient algorithms for the joint 
design of the transmit oovarianoe matrix and the opportunistio user seleetion polioy in multioell 
networks with multiple multieast groups. Here, group-rate (whieh is defined as the physioal layer 
oode-rate in a oertain bloek multiplied by the fraetion of users that ean suooessfully deeode) is 
utilized as the optimizing oriterion. This is different from most works in the literature on physioal 
layer multioasting, e.g., [4], [8], whieh typioally do not eonsider user seleetion and, thus, ean 
utilize the signal-to-interferenee-plus-noise ratio (SINK) as the optimizing eriterion. In faot, the 
SINK oriterion is not applioable when user seleetion is oonsidered sinoe it does not refleot the 
effeet of the number of users served in eaeh bloek. For example, if the system is optimized by 
maximizing the worst SINK among all seleoted users, then only one user (i.e., the user with the 
best effeotive ohannel) in eaeh group would be seleoted. However, finding the optimal transmit 
oovarianee matrix and the opportunistio user seleetion (OUS) poliey to maximize the group-rate 
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may be difficult and is, in fact, claimed to be NP-hard in [22]. Therefore, we propose in this 
work an approximate solution based on the introduction and relaxation of a set of binary user 
selection variables. This technique is similar to that previously adopted in [24] and [25] for 
admission control problems. 

In this work, we consider both single-group and multi-group multicasting scenarios. We show 
that OUS can be especially effective in the latter case, where it not only can help avoid limitations 
by the worst user but also can help coordinate interference among different cells and multicast 
groups. In the single group scenario, the optimization of the transmit covariance matrix and the 
user subset selection is performed first by relaxing the integer constraints corresponding to the 
user selection variables and then by performing a sequential deflation technique, where users are 
eliminated one-by-one from an initial set of all users to yield candidate user subsets. In the multi¬ 
group scenario, the design parameters are determined based on the group-rate balancing criterion 
(also referred to as the max-min weighted group-rate criterion) and the optimization problem 
is solved using a similar relaxation technique along with a successive convex approximation 
(SCA) approach. In the above, the channel state information (CSI) of all users is first assumed 
to be available at the transmitter, which may be costly in practice. To reduce the transmission 
overhead, we also consider the case where only a subset of users feeds back their CSI in each 
block and propose a design based on the balancing of the expected group-rates. In addition to 
the SCA approach mentioned above, a sample average approximation (SAA) technique [26], 
[27] is further introduced to handle the probabilistic terms that may arise. The effectiveness of 
the proposed schemes, compared to [21]-[23], is demonstrated through computer simulations. 

In multicell systems, three cases are often considered based on different levels of BS cooper¬ 
ation [28]-[30], namely, full BS cooperation, interference coordination, and no BS cooperation. 
Full BS cooperation refers to the case where all BSs have knowledge of the information intended 
to all multicast groups and transmit cooperatively as a networked multiple-input multiple-output 
(MIMO) system. Interference coordination refers to the case where each BS serves only one 
multicast group and only has knowledge of the message intended for that group. In this case, 
cooperative transmission of common data is not possible, but the transmit signals can be designed 
to reduce cochannel interference among different multicast groups. No BS cooperation refers to 
the case where no CSI from users in other cells is available and thus no cooperation is adopted. 
Our system model is general enough to include all the cases mentioned above. Moreover, we 
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Fig. 2. Illustration of a multicell network with 3 multicast groups. 

would also like to mention that, in this work, we do not restriet ourselves to rank-1 transmissions, 
as done in [4], [5], [8]-[10], [22], [23], in order to maintain optimality [7] and also to avoid 
distraeting the readers, sinee a rank-1 eonstraint on the transmit eovarianee matrix is already 
enough to make the problem intractable [4], even without considering user selection. However, if 
rank-1 transmissions are desired, our solution can be used as the semi-definite relaxation (SDR) 
of the optimal beamformer design and the desired beamforming vectors can be extracted from 
our solutions using rank-1 approximation techniques, e.g., the Gaussian randomization procedure, 
as discussed in [31]. 

The remainder of the paper is organized as follows. In Section II, we introduce the general 
multicell multicast scenario and the proposed group-rate balancing problem. In Sections III and 
IV, the joint design of the transmit covariance matrix and the user selection policy is examined 
for the single-group and the multi-group scenarios, respectively. In Section V, the problem is 
extended to the case where only a subset of users feeds back their CSI in each block. Finally, 
computer simulations are provided in Section VI and the paper is concluded in Section VII. 

II. System Model and Problem Formulation 

Let us consider a downlink multicell network with i? BSs, each equipped with M antennas, 
and K single-antenna users. An illustration of the system under consideration is given in Fig. 2. 
The set of users /C, with cardinality |/C| = K, is divided into G multicast groups, denoted 
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by /Cl, .. K,g, where users in the same group are interested in reeeiving the same multieast 
message. Moreover, let Eg {1,... ,B} be the set of BSs that are serving group Kg and let 
Gb ^ {1) • • •) G} be the set of multieast groups served by BS b. Notiee that the ease of full BS 
eooperation ean be considered by setting Eg = {1,, B}, for all g, and Gb = {1, • • •, G}, for 
all b, whereas the case of interference coordination, where each BS has information intended for 
only one multicast group, can be considered by setting B = G, Eg = {g} and Gh = {&}, for all 
g and b. The case of no cooperation can also be considered by setting B = G = 1, Ei = {1} 
and Gi = {!}, where the BS focuses on sending information to the intended users without 
knowledge of other BSs’ channels or transmitted signals (i.e., while neglecting the presence of 
other BSs). 

Here, we consider a delay-tolerant application where multicast messages can be encoded over 
multiple coherence intervals (hereafter referred to as blocks). Let Sfe g[n] be the signal transmitted 
by BS 6 e Eg to users in group g in block n. In this case, the received signal at user k e Kg 
can be written as 


yk[n\ = ^ hb,k[n\^Sb,g[n\ + ^ hb,k[n]^ ^ Sb/[n] + ^ hb^,k[n]^ ^ Sb'/[n] + nk[n], ( 1 ) 

b&Bg bGBg eeGt, V^Bg 

where hfe fc[n] is the M x 1 complex channel vector between BS b and user k in block n and 
nk[n] is the complex additive white Gaussian noise (AWGN) at user k with zero mean and 
unit variance, i.e., Uk ~ CJ\f{0, 1). The entries of ^ are assumed to be independent complex 
Gaussian random variables with variances that may be different for different b and k, due to 
path loss or other large-scale fading effects. The first term in (1) represents the sum of signals 
intended for user k, the second term represents the interference caused by the signals transmitted 
by user /c’s serving BSs to users in other multicast groups, and the third term represents the 
interference caused by non-serving BSs. 

By defining hk[n] = [hi,fc[n]^,..., and = [sGg[n]^,... ,SB,g[n]^]^, where 

= 0 for b ^ Eg, the received signal at user k can be rewritten as 

yk[n] = hfc [n]^Sg[n] + [n ]^ ^ [n] + nk[n]. (2) 

Therefore, the SINK at user G in block n is given by 


SINR^n] 


E[|hfc[?7,]^Sg[?7,]p] 
E£^gE[|h4n]^s4n]|2] + 1 


tr{Qg[n]hk[n]hk[n]^) 


( 3 ) 
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where Qg[ri] = E[sg[?7,]sg[?7,]^]. Note that Qg[n] is an MB x MB matrix whose (6, &')-th bloek 
entry of dimension M x M is {Qg[n]}bfi/ = E[sf,g[n]s(,/g[n]^]. This bloek entry is zero if either 
b or b' does not belong to Bg. The signals transmitted by BS b must satisfy the per BS power 
eonstraint 


^ Pb- (4) 

g^Gb g&Gb 


By employing OUS, only a subset of users in eaeh group, i.e., Ag[n\ C ]Cg, for group g, is 
required to deeode suceessfully in bloek n. In this ease, the rate of the physieal layer eneoding 
(i.e., the inner eode) of group g, i.e., Rg[n\, only needs to be low enough to ensure sueeessful 
decoding by users in Ag[n]. More specifically, by assuming that the duration of each block is 
large enough to invoke Shannon’s random coding argument [32], the rate Rg[n] of the message 


intended for the subset of users Ag[n] in block n must be chosen such that 


Rg [n] < log (1 + SINRfc [n ]) = log 14 


tr(Qg[n]hfc[n]hfc[n 


H\ 


Y^i^g tr(Qf [n ]hfc [n][n]^) +1 


(bps/Hz) (5) 


for all k E Ag[n]. 

Since the codeword in each block is decoded only by a subset of users, each user effectively 
experiences an erasure channel across multiple blocks and, thus, an outer (erasure) code, such as 
LT, Raptor, and Fountain codes [19], [20], is needed to ensure eventual recovery of the original 
multicast data at all users. By assuming that an ideal erasure code is employed so that perfect 
erasure correction is performed, the average rate of user k E fCg over T blocks can be written 
as 



n=l 


where Ip} is the indicator function, and the achievable multicast rate of group g is thus given 
by 


R 


p,min 


mm — 

fce/Cg T 


T 

^ ^ Rg['bb\'^{k£Agln]}- 

n=l 


(7) 


Conventional physical layer multicast problems often use Rg^min as the maximization criterion 
for system design. However, maximizing Rg^rnin would require non-causal knowledge of the 
channel realizations in the future T fading blocks, which is not attainable in practice. In order 
to implement the user selection and transmit signal design in real-time, we propose to maximize 
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the (average) group-rate defined as 



(8) 


for group g, where 


:Rg[n\ ^{k&Ag[n\} 
' 3' halC ^ 


(9) 


is the instantaneous group-rate in block n. The average group-rate provides a measure of the 
average speed for which the users are able to acquire the multicast message. One can see that 
the maximization of the average group-rate can be achieved by maximizing the instantaneous 
group-rate block-by-block, and hence can be implemented in real-time. However, it should be 
noted that maximizing the average group-rate Rg may not be equivalent to maximizing the actual 
multicast rate Rg,min in the general case. Therefore, Rg should only be viewed as an auxiliary 
maximization criterion that is used to allow for real-time optimization of the system parameters. 
However, it was shown in [22] (and observed similarly for single antenna scenarios in [15] and 
[18]) that, when the channel vectors are independent and identically distributed (i.i.d.) across 
users and over time, the multicast rate Rg,min is asymptotically equal to the average group-rate 
Rg and, thus, is maximized asymptotically by maximizing the instantaneous group-rate Rg[n] 
in each block. The i.i.d. assumption is reasonable in many cases and is often considered when 
prior knowledge of the channel distributions or the users’ locations are unavailable. When the 
channel vectors are non-i.i.d., normalization of the channel coefficients can be performed to 
ensure long-term fairness, as to be discussed in Sections III-B and IV-B. 

Maximizing the instantaneous group-rate for block n involves optimizing over the transmit 
covariance matrices {Qg[n]}^=i, the user subsets {Ag[n]\^^^, and the physical layer code¬ 
rates {Rg[n]]^=i subject to the power and rate constraints given in (4) and (5). This problem 
is examined for both the single-group and the multi-group multicasting scenarios. Since the 
optimization is performed separately in each block, we shall omit the block index n in the 
remainder of this work. 

In the single-group multicasting scenario (i.e., the case where G = 1), the optimal system 
parameters can be determined by solving the following instantaneous group-rate maximization 
(GRM) problem. 
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Group-Rate Maximization (GRM) Problem: 


maximize 



(10a) 

subject to 

log 2 [l +tr(Qhfchf)] 

> P, for k e A, 

(10b) 


tr({Q}M) < Pb, V6, 

Q A 0. 

(10c) 

variables: 





Notice that no cochannel interference exists in this case and that the group index g is omitted 
since only a single multicast group is considered. This problem is claimed to be NP-hard in [22] 
and no known solutions are available to solve this problem efficiently. An approximate solution 
will be studied in Section III based on the introduction of a set of binary user selection variables 
and a convex relaxation of the problem. 

In the multigroup multicasting scenario (i.e., the case where G > 1), we consider the following 
group-rate balancing (GRB) optimization problem where the system parameters are jointly 
determined to maximize the worst weighted group-rate among all multicast groups. 

Group-Rate Balancing (GRB) Problem: 


maximize 


(11a) 

subject to 

" E JwV 0 ^ 

(11b) 


G 

X^tr({Qg}fe,6) < Pb, V6, Qg A 0, 

9=1 

(11c) 


{Qg}b,b' = 0MxM, for Bg or b' ^ Eg, 

(lid) 

variables: 

{Q3}^=1, {Ag}^^l, {Rg}^^l, 



where Tg is the parameter that specifies the priority of group ICg. Solving the above problem 
guarantees a common level of quality-of-service (QoS) for all multicast groups. In particular, if 
the resulting objective value is a, then the group-rate of group g will be at least r^a, for all 
g. It is worthwhile to note that the GRB problem is analogous to the SINR-balancing problem 
often considered for the design of multiuser downlink beamforming [33] or multigroup multicast 
beamforming [8] schemes. While the SINR criterion may be suitable for conventional multicast 
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beamforming designs, where all users are required to deeode successfully in each fading block, 
it is not suitable for systems employing OUS. This is because, when considering the SINR- 
balancing (or, equivalently, the max-min SINK) criterion, the optimization will result in a trivial 
user selection policy where only the user with the best channel in each group is chosen. Even 
though the physical layer code-rate can be chosen to be the highest in this case, the group-rate 
is limited by the fact that only one user is able to decode. The GRB formulation is interesting 
in the sense that it allows us to take into consideration the effect of user selection as well as 
to incorporate the impact of QoS requirements, i.e., Tg, into the problem. This problem will 
be examined further in Section IV and an approximate solution will be proposed based on a 
relaxation similar to that in the single-group scenario. 

Remark 1: Most works in the literature on physical layer multicasting (e.g., [4]-[6], [8]- 
[10], [13], [14]) assume that all users must decode successfully in each block. In this case, the 
outer erasure code is not needed and the instantaneous group-rate will be equal to the physical 
layer code-rate in each block. The group-rate in block n must then satisfy 

Rg[n] = Rgln] < log (1 + SINR^n]), Vfc G ICg. 

The conventional approach can be viewed as the case where we choose Ag[n] = Kg, for all 
g and n. The proposed OUS scheme instead optimizes over Ag[n] and, thus, the achievable 
group-rate can be no less than that in conventional systems. 

III. Group-Rate Maximization for Single-Group Multicasting 

In this section, we examine the GRM problem in (10) for the single-group multicasting 
scenario. Note that the joint design of all system parameters in this problem is in general NP- 
hard [22] due to the combinatorial nature of the user subset selection. However, as similarly 
observed in [22], the problem reduces to a standard SDP problem when the user subset is given. 
In particular, for a given user subset A, the SDP problem can be formulated as follows: 


max 

R-\A\ 

(12a) 

subject to 

logail -f tr(Qhfchf)] > R, for k e A, 

(12b) 


tr({Q}b,b) < Pb, V6, Q ^ 0, 

(12c) 

variables: 

Q,R. 
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Due to this observation, the authors in [22] proposed a heuristic subset search algorithm where 
all users are considered initially and one user is excluded from the set in each iteration until no 
further gain in the objective value is obtained. The user removed in each iteration is the user 
whose removal from the set results in the maximum increase in the objective value. Let g{A) be 
the optimal objective value of (12) for a given user subset A. Then, the subset search algorithm 
can be summarized as follows. 

Algorithm 1 Subset Search Algorithm [22] 

1) Set ^ JL}. 

2) Define A-k ^ A — {k}, for all k E A, and let k* = argmax^ g{A-k)- 

3) If^7(^) > g{A-k*), then stop and take A as the solution; else, set AA — {k*} and go 
to Step 2. 


Notice from Algorithm 1 that, in evaluating k* in each iteration, the SDP problem in (12) 
needs to be solved for every possible choice of A-k- This requires one to solve the above SDP 
problem in the order of 0{K'^) times, which can be inefficient as K increases. 


A. Reformulation and Relaxation of the GRM Problem using User Selection Variables 


To solve the GRM problem in (10) more efficiently, let us first introduce a set of binary 
user selection variables {sfc,V/c E /C}, where = 1 if user k is selected (i.e., if k E A) and 
Sk = 0, otherwise. By choosing 5 to be sufficiently small, the GRM problem can be equivalently 
formulated as 


max 


' ' k&K 


subject to log 2 [l + tr(Qhfchf )] + (5“^(1 — su) > R, Vk E /C, 
tr({Q}fe,fe) < Pb, V6, Q ^ 0, 

Sk E {0,1}, V/c, 
variables: Q , R, {sfc}fcg/c. 


(13a) 

(13b) 

(13c) 

(13d) 


The equivalence of the two problems is stated in the following lemma. The proof is given in 
Appendix A. 
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Lemma 1: For 5 < 


maxfc log2(l + YJb=i -PfellK 


-1 


k}k&K) is an optimal 


solution of the problem in (13) if and only if {Q*, R*, A*), where A* = {k E K, ■. s\ = 1}, is 
an optimal solution of the GRM problem in (10). 


Notice that the problem in (13) is non-convex due to the integer constraints on s^’s in (13d). To 
obtain an efficient solution, we consider a relaxation of the problem where the integer eonstraints 
on Sfc’s are replaeed with the linear eonstraints 0 < < 1, for all k. By taking the logarithm 

of the objeetive and by omitting the irrelevant variables, the relaxed problem can be written as 


max lni? + ln^^Sfc ( 14 a) 

keic 

subject to log 2 [l + tr(Qhfch^)] + (5“^(1 — Sfc) > i?, \/k e K., (14b) 

tr({Q}fe,fe) < Pfe, V6, Q P 0, (14c) 

0<Sfc<l,VA:, (14d) 

variables: Q , R, {Sk}kelC- 


This problem is eonvex and ean be solved efficiently using general purpose solvers sueh as CVX 
[34]. By relaxing the integer constraints on the user seleetion variables, the term (5“^(1 — Sk) 
can take on any value between 0 and 5“^, and ean be viewed as a measure of rate violation, i.e., 
the difference between the code-rate R and the achievable rate of user k. Similar approaches 
have also been used in studies of admission control problems in [24] and [25]. Notice that the 
solution of Sk in the relaxed problem may not take on the value 0 or 1. To eonvert the solution of 
the relaxed problem to a feasible solution of (13), we propose a sequential deflation teehnique, 
where the D users with the D smallest values of Sk are removed in each iteration, after which 
the values of Q, R, and Sk, Vfc, are updated by solving (14) again. The value of D can be ehosen 
as 1 in most cases but ean be chosen to be greater than one for eomplexity reduetion. Among 
the \K/D~\ possible subsets obtained from this procedure, where [a] is the smallest integer not 
less than a, the user subset that yields the largest group-rate is ehosen. The sequential deflation 
algorithm is summarized in Algorithm 2. 

Notiee that the sequential deflation algorithm proposed above is different from that proposed 
in [25] where the algorithm is to terminate whenever the removal of a user no longer results in 
an increase in group-rate. However, when applying their approach to our problem, the algorithm 
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Algorithm 2 Opportunistic User Selection by Sequential Deflation 

(i) Initialize by setting ^ /C, 0, and a* <(— 0. 

(ii) Solve (12) for given user subset A and let a be the resulting objective value. 

(iii) If d > a*, then set a* d and A* A. 

(iv) Solve the relaxed problem (14) for /C = ^ to yield the values of Sk for all k e A. 

(v) Set ^ ^ ^ — iSmin, where Smin is the set of users associated with the D smallest values 
of Sk among users in A. 

(vi) Repeat steps (ii)-(v) until A = 0. Then, take A* as the desired user subset. 


often terminates early in the process at a local optimum that is close to choosing all users 
as the serving subset. Moreover, it is worthwhile to remark that, in Algorithm 2, two convex 
optimization problems are solved in each iteration (i.e., one for computing the relaxed problem 
(14) and one for solving (12) for given A) and the number of iterations is equal to the number 
of users K in the worst case, when D = 1. This is a significant improvement over the 0{K‘^) 
worst-case complexity required in the subset search method summarized in Algorithm 1 [22]. 

B. Fairness of the GRM Problem in the Non-I.I.D. Case based on Channel Normalization 

Notice that maximizing the instantaneous group-rate, as done in the previous subsection, 
equivalently maximizes the average group-rate in (8). When the channel vectors (i.e., the short¬ 
term fading coefficients) are i.i.d. across users, all users will have equal probability of being 
selected in each block and, thus, the average group-rate would eventually be the rate achieved 
by all users as T ^ cxo [22]. However, when the channel vectors are non-identically distributed, 
users with better average channels (namely, users closer to the BSs) may have a higher probability 
of being selected under our proposed scheme. In this case, an issue of fairness may arise, similar 
to that observed in [15] and [18] for the single-antenna case. In this section, a heuristic channel 
normalization technique is proposed to address the aforementioned fairness issue. 

The key idea is to normalize the channel vectors of the users by their long term statistics and 
perform the proposed OUS based on the normalized channel vectors (which will then be i.i.d.). 
More specifically, suppose that the channel vector hb^k[n] between BS b and user k in block n 
has entries that are i.i.d. CAf{0, d^^), where db,k is the distance between BS b and user k and a 
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is the path loss exponent. Moreover, let (ife = ZI/cga: ^b,k be the average distanee of all users 
to BS b. Let us define the normalized ehannel veetor between BS b and user k in bloek n as 



(15) 


whieh has entries that are i.i.d. CAf{0, “). The normalized ehannel veetors are then utilized to 


eompute the optimal user subset, denoted by A*, using the algorithm proposed in the previous 


subseetion. The optimal transmit eovarianee matrix ean then be computed by solving (12) for 


given A* using the original channel vectors hf, fc, V6, k. 


The normalized channel vectors are scaled by their average distance so that the distance 
between the BSs and the users are taken into account in the signal-to-noise ratio (SNR) at the 
receiver. After normalization, the channel vectors experienced by all users will become i.i.d. 
and, by performing OUS based on the normalized channel vectors, all users will have equal 
probability of being selected in each block, ensuring fairness among users. The users that are 
opportunistically selected in each block will likely be users whose instantaneous channel gains 
are larger than their respective averages, exploiting the advantages of both temporal and multiuser 
diversity. It is also interesting to note that, even though the users in the subset A* are chosen 
as target users and the transmit covariance matrix Q is chosen to maximize the rate of users 
in A*, it is possible that users outside of the target subset A* may also be served due to their 
large average channel gains (i.e., their close distance to the BSs) or because of their coincidental 
locations in the directions of the beamformed signals. The performance of the proposed fair 
OUS scheme is demonstrated through simulations in Section VI. 

It is necessary to note that the proposed normalization scheme is applicable to any case in 
which the average channel gain is different for different users, not limited to that caused by path 
loss. For frequency-selective fading scenarios, multicarrier modulation schemes such as OFDM 
can be considered and the proposed algorithm, including the channel normalization technique, 
can be applied to each subcarrier individually. 


IV. Group-Rate Balancing for Multigroup Multicasting 


In this section, we examine in more detail the GRB problem described in (11) for the 
multigroup multicasting scenario. The problem is more challenging than that in the single-group 
scenario because of the existence of cochannel interference among different groups. 
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Notice that, similar to the single-group scenario, the complexity of solving the GRB problem 
in ( 11 ) is largely due to the eombinatorial nature of the search for the optimal user subsets 
{'^g}f=i- However, different from the problem in (12), the optimization problem does not reduce 
to a eonvex optimization problem even when the user subsets are given. Specifieally, 

when are given, we have 

max min — R„ (16a) 

g&{l,...,G} T'g 

subject to log, (l + + ^0' 

G 

X]tr({Q 3 }fc, 6 ) < Pb, V 6 , Qg h 0, yg, (16c) 

9=1 

{Qg}fe,6' =0 mxm, for b^Bg or b' ^ Bg, (16d) 

variables: 

where r' = Tg\}Cg\/\Ag\. We ean observe that the rate variables {Rg}^^i are uneonstrained 
below and that the objective value does not become smaller by decreasing Rg/r'g to the value 
m\ng^{i^,„^G}Rg/'Tg for AX g = 1,... ,G. Consequently, we ean impose the eonstraint 

Pi/'T'i = ■■■ = Rg/t’g 

on problem (16) without loss of optimality. 

By introducing the variable a = Ri/t[ = ■■■ = Rg/t'q, the problem in (16) can be 
equivalently expressed as 

max a (17a) 

subjeet to tr(Qghfchf) > j ^tr(Q£hfchf) + 1 j , V/c e Ag^^g, (17b) 

^ \i-+g ) 

G 

^tr({Qg}b,b) < Pft, V6, QgP0,V5(, (17o) 

9=1 

{Qg}b,b' = 0MxM, for b^Bg orb' ^ Bg, (17d) 

variables: {Qg}^=i, a. 
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This problem is still non-convex due to the eonstraint in (17b). However, for a fixed a, the 
eonstraint beeomes eonvex and the problem beeomes a eonvex feasibility problem. The optimal 
value of a ean then be found via a biseetion seareh as deseribed in the following. 

To perform the biseetion seareh, it is neeessary to first determine an upper bound and a lower 
bound on the value of a. To do so, notiee that, in (17), r'^a = Tga\lCg\/\Ag\ ean be viewed as 
the minimum rate aehievable by all users in the subset Ag and, thus, is upper-bounded by the 
rate aehievable when all BS powers are used to beamform data to the best user in group g (and 
no signal is sent to other groups), i.e.. 



(18) 


where ^ is the ehannel from BS b to user k. This inequality must hold for all g and, thus, a 
is upper-bounded by 



(19) 


Moreover, sinee the rate is non-negative, a is trivially lower-bounded by 0. The biseetion 
algorithm repeatedly biseets the interval between the upper and lower bounds given above until 
the solution is obtained. Details are given in Algorithm 3. 

Algorithm 3 Biseetion Seareh for a 

1) Initialize (7 as in (19) and L = 0. 

2) Set and solve the feasibility problem eorresponding to (17). If it is feasible, then 

set L = a. Otherwise, set U = a. 

3) Repeat Step 2 until U — L < e, where e is the value speeifying the eonvergenee eriteria. 
The resulting {Qg}^=i and a are the desired solution to (17). 


In the above, the optimal transmit eovarianee matriees {Qg}^=i (and the eorresponding rates 
{Rg}^=i) were found for given user subsets {Ag}^^^. To find {Ag}^^i, one ean employ a 
heuristie OUS poliey similar to the subset seareh algorithm deseribed in Algorithm 1. However, 
as mentioned previously, this requires solving (17) in the order of O(iT^) times, with a biseetion 
seareh embedded in eaeh eomputation. In the following, we propose a more effieient method 
based on the introduetion of binary user seleetion variables, as done in the single-group seenario. 
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Here, an SCA approaeh is further adopted to eope with the non-eonvexity eaused by the 
interferenee terms. 


A. Reformulation and Relaxation of the GRB Problem using User Selection Variables 


Speeifieally, following the teehnique given in Seetion III, let us introduee the set of binary user 
seleetion variables {sfc, V/Cg, V^f}, where = 1 if /c G for some g, and Sk = 0, otherwise. 
Then, for a suffieiently small 5, the GRB problem in (11) ean be equivalently formulated as 


max mm 


ge{l,...,G} Tg\lCg 


kelCn 


Sk 


subjeet to log2 1 + 


tr(Qghfchf) 


tr(Q£hfch|^) + l 


+ 6 (1-Sfc) > Rg, Vk e ICg, Vg, 


G 


< Pfe, V6, Qg t 0,Vg, 

9=1 


{Qg}b,b' = 0MxM, for h ^ Eg or U ^ Eg, 
Sk G {0,1}, V/c G /C, 
variables: {Qg}^=i, {Rg}g=i, 


(20a) 

(20b) 

(20e) 

(20d) 

(20e) 


The equivalenee of the problems in (11) and (20) is stated in the following lemma. 
Lemma 2: For 


5 < 


B 

max max logo ( 1 + > Pb\\hbk 

g&g keiCg II , 

0=1 


-1 


{Sk}keK.) is an optimal solution of the problem in (20) if and only 
{R*g}^^i, {A*g}^^i), where A*g ^ {k E ICg : si = 1}, for g = 1,... ,G, is an optimal solution of 
the GRB problem in (11). 


The proof is similar to that of Lemma 1 and, thus, is omitted. Moreover, let us also eonsider 
a relaxation where the integer eonstraints on sRs are replaeed with the linear eonstraints 0 < 
Sk < for all k. By introdueing the auxiliary variable a, the relaxed problem ean be written in 
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the epigraph form 


max a. 

(21a) 

subject to Rg^ Sk> rg|/Cg|a^ ^g, 

(21b) 

k&Kg 


^fc({Q£}f=i) + - Sk) > Rg,yk e /Cg, 

(21c) 

0 < Sfc < 1, Vfc G /C, (20c), and (20d), 

(21d) 


variables: {Rg}%i, {sfc}fcG^, «, 


where 

rk{{Qe}f=i) = logs 

for k E K,g. Notice that this problem cannot be solved using the bisection search algorithm 
since, even when a is given, the problem does not reduce to a convex feasibility problem. 
This is because the value of Rg cannot be determined explicitly and the constraint in (21c) is 
non-convex. To address this issue, we employ an SCA approach as described in the following. 
First, notice that, since Rg is positive, the non-convex constraint in (21b) can be written as 

Sk - a\^lTg\lCg\R~^a\^lTg\lCg\ > 0. (23) 

keKg 

Then, by applying the Schur complement [35], this constraint can be equivalently written as the 
linear matrix inequality constraint 

^^ ^ 0. (24) 

Secondly, to address the non-convexity of the constraint in (21c), we adopt an SCA technique 
similar to that employed in [36]. Specifically, given any satisfying (20c) and (20d), 

which together with a = Rg = 5“^, Sk = 0, k, is a feasible point to (21), we can rewrite 


1 + 


tr(Qgh, 


.hf) 


tr(Q£hfch|^) + 1 


( 22 ) 
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and obtain its lower bound as 

rk{{Qi}^=i) 


= logs 1 

l + ^tr(Qfhfchf)] 

-loga 1 


\ e=i J 

V 


( ^ ^ 

, / 

> log21 

1 +^tr(Q£hfchf) 
\ i=\ ) 

) -iog2 h 






if) 


if) 


rkme}?=i I {Qaf=i) 


J2i^g tr((Q£ - Q£)hfch|^) 
[1 + E£^3tr(Qfhfchf)]ln2 


where the inequality eomes from the first-order eondition of the eoneave funetion log 2 (-)> i-e., 
logal/ < \ 0 g 2 X + - x) for any x,y > Q. Note that fk{{Qi}f^^ \ is eoneave 

in By replaeing rk{{Q,i}^=i) in eonstraint (21e) with \ and by 

replaeing (21b) with (24), the problem in (21) ean be approximated as 


max a (25a) 

subjeet to ( ^\ ^25b) 

V EfcGX.Sfc / 

rk{{Q,i}^=i I {Q^}f=i) + 5~^0--Sk)>Rg, Vfc G /Cg, Vfif, (25o) 

0 < Sfc < 1, VA: e /C, (20o), and (20d), (25d) 

variables: {Qg}f=i, {Rg}%i, {sfc}fcGX, «• 

The resulting optimization problem in (25) is eonvex and ean be solved using general purpose 
solvers sueh as CVX [34]. Notiee that this approximation is eonservative in the sense that solving 
(25) yields a feasible approximate solution to (21). However, the approximation performanee ean 
be improved by iteratively applying the same approximation to (21), with taken as the 

solution of obtained in the previous iteration. The proeess ean be repeated until the 

value of OL eonverges. The algorithm is summarized in Algorithm 4. 


Algorithm 4 Relaxed GRB via Sueeessive Convex Approximation (SCA) 

1) Initialize by any feasible point of (20e), (20d), and set (5 = 0. 

2) Solve problem (25), and denote the solution of as and the objeetive value as a*. 

3) If \a — a*\ < e, then stop; else, set •<— for all i, a ^ a*, and go to Step 2. 
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The convergence of Algorithm 4 is described in the following proposition. 

Proposition 1: Let {{Qg}^=i, denote the solution of (25). Then, 

any limit point of the iterates generated by Algorithm 4 is 

a stationary point of problem (21). 

The proof is given in Appendix B. After obtaining the values of from the relaxed problem 
(and the SC A approach in Algorithm 4), the sequential deflation technique proposed in the 
previous section can then be utilized to obtain a solution for the user subsets {Ag}'f^i. Similarly, 
by employing the sequential deflation technique, the optimization problem in (25) is solved 0{K) 
times, which is a significant improvement over the subset search with bisection. The performance 
of the two algorithms will be compared in Section VI. 

B. Fairness of the GRB Problem in the Non-l.l.D. Case based on Channel Normalization 

Notice that, similar to the single-group scenario, the GRB problem may also suffer fairness 
issues since the average group-rate is again used as the maximization criterion in each group. 
In this section, we extend the heuristic channel normalization technique, previously proposed in 
Section III-B, to the multi-group scenario. 

Specifically, suppose that the channel vector lib fc[n] between BS b and user k in block n has 
entries that are i.i.d. CAA(0,(i^^), where db^k is the distance between BS b and user k and a is 
the path loss exponent. Let db,g = J^k&iCg be the average distance between BS b and 
users in group g. Then, the normalized channel vector between BS b and user k in group g can 
be defined as 



(26) 


Here, instead of scaling the normalized channel vectors by the average distance associated with 
all users, as in (26), the channel vectors are scaled by the average distance of users within 
its own group. Similar to Section III-B, the normalized channel vectors are then utilized to 
compute the optimal user subsets, denoted by A*g, for all g, using the algorithm proposed in 


the previous subsection. With user subsets {A*}^^i, the optimal transmit covariance matrix can 


then be computed by solving (16) using the original channel vectors fc, V6, k. 

By performing OUS based on the normalized channel vectors, users in the same group will 
have equal probability of being selected in each block. These users will likely be those whose 
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instantaneous channels are temporarily more favorable, either in the sense that their channel 
gains are large relative to their respective averages or in the sense that their channel directions 
cause less interference to other groups. The performance of the proposed fair OUS scheme is 
demonstrated through simulations in Section VI. 

V. Group-Rate Balancing based on Partial Channel Feedback 

In previous sections, the design of the transmit covariance matrices and the user subsets was 
performed based on instantaneous knowledge of all users’ CSI. In practice, this may require a 
large amount of feedback, especially when the number of users in the system is large. In this 
section, we consider the case where only a subset of users feeds back their CSI in each block. 
Since the CSI of some users may be missing, the optimization can only be performed based 
on the expected group-rate, i.e., the expectation of the average group-rate with respect to the 
statistics of the unknown channels. 

Let /Cg_fb be the set of users from group K,g that have chosen to feedback their channels in 
the given block and let = ICg — ICg^tb be the complement set, i.e., the set of users from 
group JCg whose CSI is unknown. In this case, the system parameters can be found by solving 
the group-rate balancing with partial feedback (GRB-PF) problem as given below: 

GRB with Partial Feedback (GRB-PF) Problem: 

^9 TP 

max mm — rrr-rtj 

G 

subject to ^tr({Qg}b,b) < Pt, V6, Qg P 0,^g, (27b) 

9=1 

{Qg}b,b' = OMxM, for b^Bg or b' ^ Eg, (27c) 

variables: {Qg}^^^, {Rg}^^^, 

Here, we assume that a user’s decision to feedback its channel (or not) is independent of its 
channel realizations. Otherwise, the expectation must be evaluated by conditioning on the specific 


^{r-fc({Qati)>R4 


k£Kn 


(27a) 
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feedback strategy^ In particular, the expectation term inside the objective function can be written 
as 


El 

fb 


log2 1 + 


tT(Qgh^hjf) 


E£7^gtr(q£hj.h") + 1 


V Prolog (l I \ 


— ^9 


■ (28) 


Notice that the probability in the second term cannot be evaluated in closed-form and, thus, 
is difficult to handle from an optimization perspective. To address this issue, we propose to 
consider an approximation technique where the probability is replaced with its sample average 
approximation (SAA). Similar techniques have also been widely adopted in the area of stochastic 
optimization, e.g., [26], [27]. 

More specifically, let be the density function of and let be J vectors randomly 

generated according to /h^,. Then, the probability can be approximated as 


Prlog2 1 + 


tr(Qghfch 




— ^9 


J 


El 

2=1 


log2 1-1- 


tr(Qgh» (h» )g) 


(29) 


With knowledge of the distributions of the channel vectors h^, for each k G JCg^, the objective 
function in (27) can be approximated as 

J 


Rn 


mm —— 

9 Eg I /Cg 


Ef 


1 


{r-;o({QAf=i)>«g} 
fcGXg^fb keiC' 


E 7 Ei 


J . 

g,fb i ^ 


{^■fc,j({Q£}£Ll)>12g} 


(30) 


where 

rk,ji{Qe}f=i) = log2 

for k e /C^ fb- The channel vectors can be viewed as the channel vectors of J virtual 

users associated with user k. The entire set of virtual users in group g is defined as /Cg^vir — 
/C^ fb X {1,..., J}. Moreover, let ^g,fb ^ ^g,ih and ^g,vir ^ be the subsets of selected 

feedback and virtual users, respectively, in group g that are intended to successfully decode, and 


1 + 


tr(QghSf^(hSfV) 




n + i. 


(31) 


'Note that, in cases with partial feedback (i.e., when only part of the users feedback their channels), the specific feedback 
strategy may have a significant impact on the group-rate performance. In particular, allowing channel feedback from users with 
bad channels may improve service towards these users in each block but may limit the achievable multiuser and temporal 
diversity gains over time; and allowing feedback from users with good channels may have the opposite effect. The optimal 
feedback strategy should exploit the tradeoff between these two effects, especially in the non-i.i.d. case, but requires further 
studies that go beyond the scope of this work. Readers are referred to [37] and [38] for further discussions on this topic. 
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let the index pair (A;, j) denote the j-th virtual user assoeiated with (non-feedbaek) user k. The 
GRB-PF problem in (27) ean then be approximated as follows: 


max mm ^ ( |A,fbl + JA.virM 

g€{l,...,G} Tg\}Cg\ V J J 

(32a) 

subject to VkiiQejf^i) > Rg, VA; G ^g,fb,V5(, 

(32b) 

rk,j{{Qe}i=i) > ^9) V(A:, j) G ^g,vir,V^, 

(32c) 

(27b) and (27c), 

(32d) 


variables: {Avir}^=i- 


Notiee that the above problem is diffieult to solve due to the eombinatorial nature of the seareh 
over the user subsets and However, when and {^g_vir}^=i are 

given, the problem becomes the same as (16) with r'g = T,|/Cg|/(|^ 3 ,fb| + j|-Ag,vir|) and, thus, 
can be solved using the bisection search algorithm described in Algorithm 3. 

To determine the user selection (i.e., the choice of {Ag^fo}^=i and {^g,vir}^=i), we introduce, 
for each feedback user k G /Cg,fb, a user selection variable Sk and, for each virtual user (fc, j) G 
/Cg vir, a user selection variable tkj. Then, similar to the previous sections, the approximated 
GRB-PF problem can be reformulated as 

max min^ I ^ | (33a) 

9 yfce/Cj_fb {k,j)eK.gy^^ j 

subject to rfc({Q£}f^i) + 5"^(1 - Sk) > Rg, 'ik G /Cg,fb, 'ig, (33b) 

^fcj({Q4f=i) + - R,j) > ^3) V(A:, j) G /Cg,vir, V^f, (33c) 

Sfc G {0,1}, VA: G /C 3 ,fb, Wg, (27b), (27c), (33d) 

tk,j G {0,1}, V(A;,j) G /Cg,vir, Vfif, (33e) 

variables: {s^, VA: G /Cg,fb}J=i, {4j, V(A;, j) G /Cg,vir}J=i. 

where t'J = rJ/CJ. 

By relaxing the integer constraints and by applying properties of the Schur complement, we 
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can obtain a similar relaxed problem in its epigraph form, i.e., 


max a 

/ 

subjeet to 

V 


Rn 




\ 




fce/Cg 
iG ^ 




tk,i 


^0, Wg, 


+ <5 (1-Sfc)>i?g, V/c e/Cg,fb, V^f, 

> Rg, y{k,j) e /Cg,vir, v^, 

0<Sfc<l, Vfce/C,,fb, (27b), (27e), 

0 < tkj < 1, V(/c, j) e /Cg,vir, V^f, 


(34a) 

(34b) 

(34e) 

(34d) 

(34e) 

(34f) 


variables: {QjJLi, {sfc, V/c G /Cg,fb}Ei, {4j, V(/c, j) e /Cg,vir}Ei,a 


|G 


G 


|G 


The above problem is still non-eonvex due to the eonstraints in (34e) and (34d). Therefore, 
we adopt an SCA approaeh, similar to that in the previous seetion, where the left-hand-side 
of the inequalities in (34e) and (34d) are approximated by their eoneave lower bounds. More 
speeifieally, the problem in (34) is approximated as 


max a 

/ 

subjeet to 


Rn 


Oi\ T' 


\ 




E Sfe+TEf 




^0, yg, 


k^K^g / 

^fc({Q^}Ei I {Q^lEi) + vfc e /Cg,fb, Vfif, 

rk,j{{Qe}f=l I {Q£}f=l) +(5“^(l-4j) > Rg, V(fc, j) G /Cg,vir, yg, 
0<Sk<l,yke /Cg,fb, (27b), (27o), 

0 < 4j < 1, y{k,j) G /Cg,vir, yg, 


variables: {QjEi, {s*,, Vfc G /Cg,fb}Ei, {4j, V(fc, j) G /Cg,vir}Ei, 


a, 


(35a) 

(35b) 

(35e) 

(35d) 

(35e) 

(35f) 

G 


where {Q,e}i=i can be any point satisfying eonstraints (27b) and (27o), and rk,jW^}i=i 
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{Qe}f= is defined as 


£=1 


I {Q£}^=i) =log2(l+5Ztr(Q£hi^^(h^fc^^)^)) -log2(l+^tr(Q^h^fc^^h2^ 


Ji. 






The resulting optimization problem (35) is eonvex and ean be solved using general purpose 
solvers sueh as CVX [34]. Similarly, refined approximate solutions ean be obtained by solving 
a series of eonvex optimization problems where, in eaeh iteration, are taken as the 

solutions of {Q^}^! obtained in the previous iteration. The proeess ean be repeated until the 
objeetive value a eonverges. The algorithm is summarized in Algorithm 5. 




In 2 


Algorithm 5 Relaxed GRB-PF via SAA and SCA 

1) Randomly generate ehannel veetors V/c G V^. 

2) Initialize by any feasible point of (27b), (27e), and set d = 0. 

3) Solve problem (35), and denote the solution of as , and the objeetive value as a*. 

4) If \a — a*| < e, then stop; else, set •(— Q^, a ^ a*, and go to Step 2. 


Proposition 2: Let e /Cg,fb}^=i, V(fc,j) 

denote the solution of (35). Then, any limit point of the iterates generated by Algorithm 5 is a 
stationary point of problem (34). 

The proof of Proposition 2 is similar to that of Proposition 1 and, thus, is omitted. After 
obtaining the values of Sk and j from the relaxed problem, a sequential deflation teehnique, 
similar to that in Algorithm 2, is again applied to obtain a solution for the feedbaek and virtual 
user subsets, i.e., {'^g,vir}^=i’ respeetively. The sequential deflation algorithm 

used for the GRB-PF problem is formally deseribed in Algorithm 6. Notiee that, in this problem, 
SAA introduees a large number of virtual users whieh may inerease the eomplexity of the 
sequential deflation approaeh. To reduee eomplexity, we ehoose to eliminate more than one 
virtual user at onee in eaeh iteration. Speeifieally, if the smallest user seleetion variable is asso- 
eiated with a virtual user (i.e., if min(fcj)e^^ < minlminfegUg^g,^ ^k, 
then Dt > 1 (e.g., Dt = 5) virtual users with the Dt smallest user seleetion variables are 
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eliminated from the potential virtual user subset Agy-^r- On the other hand, if < 

min{minfc 6 U 9 ^ 3 ,fb then Ds = I feedbaek user is eliminated from the 

potential feedbaek user subset Ag^po. The removal of virtual users ean be viewed as the removal 
of ehannel realizations that are not able to jointly support a high transmission rate. 


Algorithm 6 OUS by Sequential Deflation for the GRB-PF Problem 

(i) Initialize by setting Ag^^b ^ /Cg,fb, Agy„ ^ /Cg,vir, ^ 0, -^ 3 ,vir ^ 0, for 51 = 1 ,..., G, 

and a* ^ 0. 

(ii) Solve (32) for given user subsets {Ag^ft,}^^^ and {^c,,vir}^=i, and let a be the resulting 
objeetive value. 

(iii) If d > a*, then set a* a, A*g^^ ^ -4,g,fb. and A*gy^^ ^ Agy^r, ^9- 

(iv) Solve the relaxed GRB-PF problem in (35) (using Algorithm 5) for /Cg,fb = ^g,fb and 

= Ag^ir, 'ig, to yield the values of s^, 'ik e ^g,fb, and ^{k,j) ^ Ag^^ir, 'ig. 

(v) If minfce^^ jj^Sfe < min {minfc6U9.4g,fb Sfe, niin(fcj)eUg.A 3 ,vir 4,j}, then set Ag^^^ ^ 
Ag^ih — ^min.fb, where iSmin,fb is the set of feedbaek users assoeiated with the 
Dg smallest values of Sk among users in Ag^fo\ else, if mini^k,j)<^Ag vir ^k,j < 
mm{mmk(zugAg,ihSk,mm(^k,j)&UgAg,.ir'tk,j}, then set Agy,, ^ A.vir - ‘^min.vir, where 
^min.vir is the Set of virtual users assoeiated with the Dt smallest values of tkj in Ag^vh- 

(vi) Repeat steps (ii)-(v) until = 0 for some g. Then, take {^g_fb}^=i and as 

the desired set of feedbaek and virtual users. 


It is worthwhile to note that, in the proposed seheme, the virtual users are generated at random 
and, thus, the solution obtained under eertain realizations of the virtual users’ ehannels may be 
worse than the ease without virtual users (i.e., solving the GRB problem assuming only the 
existenee of the feedbaek users). This is espeeially the ease when the number of virtual users, 
i.e., J, is small. Henee, in praetiee, one ean take the better of the two sehemes (i.e., the seheme 
eonsidering only feedbaek users’ ehannels and the proposed GRB-PF seheme) in eaeh realization 
to guarantee that no loss oeeurs due to sueh randomness. This allows us to reduee the required 
number of virtual users in praetiee. Moreover, when the users’ ehannels are i.i.d., the same set 
of virtual users ean be used for all users, whieh further reduees the eomplexity. 
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Fig. 3. For the single-group multicasting scenario with B = S BSs, M = 2 antennas per BS, and if = 30 users, we show the 
average group-rates of single-group multicasting using the proposed sequential deflation, the subset search, the broadcast, and 
Low’s algorithms. 


VI. Simulations Results 

In this section, the effectiveness of the proposed schemes is demonstrated through computer 
simulations. In the experiments, we consider a multicell network with B = 3 BSs, each equipped 
with M = 2 transmit antennas, and set the power constraints of the BSs as Pi = P 2 = P 3 = P. 
The SNR is defined as P/a'^, where = 1 is the noise variance at each receiver (as chosen 
in Section II), and the entries of the channel vectors are assumed to be i.i.d. CAf{0, 1) unless 
mentioned otherwise. The results are obtained by averaging over 600 channel realizations. 

A. Single-Group Multicasting Scenario 

First, we consider the single-group multicasting scenario with B = 3 BSs serving collabora- 
tively K = 30 users. We would like to emphasize that the scenario under consideration is different 
from having a single BS with 6 antennas since each BS here is subject to their own individual 
power constraint P. In Fig. 3, we compare the average group-rate of single-group multicasting 
using the proposed sequential deflation technique (described in Algorithm 2) with that of single¬ 
group multicasting using the subset search algorithm (described in Algorithm 1), the broadcast 
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Fig. 4. For the single-group multicasting scenario with B = S BSs, M = 2 antennas per BS, and if = 30 users, we show the 
average group-rate improvements obtained with the proposed sequential deflation algorithm and the subset search algorithm. 


scheme (where all users are served in each block), and Low’s algorithm [21] (which is based on 
a heuristic semi-orthogonal user selection algorithm). Here, Low’s algorithm is performed with 
power allocation that takes into consider the individual power constraints at different BSs. We 
ean see that, in the single group scenario, the proposed and the subset-seareh based OUS policies 
perform better than the broadcast scheme and, in fact, provide more advantages in the low SNR 
regime than in the high SNR regime. This is because, in the high SNR regime, a significant 
rate loss may be experieneed when a user is eliminated and, thus, it is preferable to serve all 
users simultaneously (as done in the broadeast seheme). In eases with OUS, the subset search 
algorithm performs slightly better than the proposed scheme, but the differenee is not significant 
and comes at the cost of much higher complexity. Low’s algorithm performs the worst among 
all scehems since the precoder is not shaped in accordanee with the individual power constraints, 
but chosen to maintain orthogonality among different signal directions [21]. In Fig. 4, we show 
the group-rate improvement of the proposed and the subset search algorithms. The group-rate 
improvement is defined as the differenee in group-rate between the compared algorithm and the 
broadcast scheme, normalized by the group-rate of the latter scheme. We ean see that, at low 
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Fig. 5. For the single-group multicasting scenario with B = S BSs, M = 2 antennas per BS, and if = 30 users, we show the 
average number of selected users obtained with the proposed sequential deflation, the subset search, and Low’s algorithms. 


SNR, the group-rate improvement ean be over 20% for both the proposed and the subset seareh 
algorithms. 

In Fig. 5, we show the average number of users that are seleeted in eaeh bloek in the proposed 
sequential deflation, the subset seareh, and Low’s algorithms. We can see that the number of 
selected users increases with SNR in all schemes. The subset search algorithm eliminates fewer 
users because it terminates whenever no further improvement is obtained after removing a user 
whereas the proposed scheme first removes users sequentially until no user remains to yield 
|/C| candidate user subsets and then chooses the solution that yields the best group-rate. Even 
though the subset search algorithm is able to choose the most appropriate user to eliminate in each 
iteration (since it performs an exhaustive search among all remaining users in each iteration), it 
may have terminated prematurely because of the existence of many locally optimum solutions. 
It is worthwhile to note that the number of users selected in each block does not reflect the 
long-term fairness of the scheme. It only shows how each scheme exploits the tradeoff between 
multiuser diversity and multicast gains. When the users’s channels are i.i.d., each user has equal 
opportunity of being selected in each block and thus, with the incorporation of the outer code. 
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Fig. 6. For the single-group multicasting scenario with B = 3 BSs, M = 2 antennas per BS, and SNR=—15dB, we show the 
average group-rate of the proposed, the subset search, and the broadcast schemes as the number of users K increases. 


the average rate aehieved by all users is asymptotieally the same. 

In Fig. 6, we show the average group-rate of the proposed, the subset search, and the broadcast 
schemes with respect to the number of users, i.e., K. The receive SNR is fixed as —15dB. Due 
to the high computational complexity of the subset search scheme, its performance is averaged 
only over 300 channel realizations for K < 125 and 30 channel realizations for K > 125. We 
can see that the OUS schemes (i.e., the proposed and the subset search algorithms) perform 
better than the broadcast scheme, especially as the number of users increases. This is because, 
in the broadcast scheme, the group-rate is limited by the worst user in the group and the channel 
conditions of the worst user will degrade continuously as the number of users increases. However, 
in the OUS schemes, the rate limitations caused by the worst users are alleviated by selecting 
users with sufficiently reliable channels in each block. For large K, the subset search algorithm 
performs better than the proposed scheme in the single-group scenario, but requires significantly 
higher computational complexity and rapidly becomes intractable as K increases. 

Next, we consider the single-group multicasting scenario with B = 3 BSs serving collab- 
oratively K = 30 users with non-identically distributed channel vectors. We assume that the 
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Fig. 7. For the single-group multicasting scenario with B = 3 BSs, M = 2 antennas per BS, and K — 30 non-i.i.d. 
users, we show the average group-rates of single-group multicasting using the proposed OUS scheme with and without fairness 
considerations, and the broadcast scheme. The minimum and 90% rates of the proposed fair OUS scheme is also shown for 
comparison. 


users are uniformly distributed in a [0, 800]mx [0, 733]m region and the eoordinates of the 3 
BSs are (150,150), (650,150), and (400,583), respeetively. The BS loeations are ehosen sueh 
that they are distaneed equally by 500 meters. The ehannel veetor hfe fc[n] is assumed to have 
entries that are i.i.d. CAf{0, with path loss exponent a = 2.5. Here, 6 sets of random user 
loeations are eonsidered, eaeh averaged over 100 ehannel realizations. In Fig. 7, we show the 
average group-rates versus the transmit SNR of the proposed single-group multieasting seheme 
with and without fairness eonsiderations as well as that of the broadeast seheme. The minimum 
rate among all users (i.e., (7)) as well as the rate aehieved by 90% of users (ealled the 90% 
rate) are also shown for eomparison. The transmit SNR is given by the transmit power P since 
the noise variance is set as 1. The proposed scheme with fairness considerations refers to the 
scheme that utilizes normalized channel vectors to compute the user subsets whereas the scheme 
without fairness refers to the scheme that utilizes the original channel vectors to compute the 
user subsets. We can see that the proposed fair OUS scheme can still achieve average group-rate 
that is higher than the broadcast scheme even though a loss is experienced compared to the 
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Fig. 8. For the single-group multicasting scenario with B — 3 BSs, M = 2 antennas per BS, and K = 3G non-i.i.d. users, 
we show the average number of selected and served users obtained with the proposed OUS scheme with and without fairness 
considerations. 


case without fairness. The minimum rate, however, ean be lower than the rate achieved by the 
broadeast seheme if the message is not transmitted over a suffieiently large number of ehannel 
realizations. By averaging over 100 ehannel realization in this figure, the minimum rate is slightly 
lower than the broadeast seheme due to the diversity of the rates among users. However, at least 
90% of users experienee rates higher than that of the broadeast scheme, even though our seheme 
is derived using the average group-rate (instead of the minimum rate in (7)) as the optimization 
eriterion. This demonstrates the effeetiveness of the proposed normalization. The minimum rate 
will beeome eloser to the average rate as the transmission oeeurs over larger number of ehannel 
realizations. 

In Fig. 8, we show the average number of users that are seleeted in the sehemes with and 
without fairness. We ean see that, in the case without fairness eonsiderations, less users are 
seleeted, and the seleeted users almost always eorrespond to users elose to the BSs. As mentioned 
in Seetion III-B, in the proposed fair OUS seheme, the users that are selected as target users are 
not the only users that may actually be served. In faet, by setting the rate to satisfy the worst 
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Fig. 9. For the multi-group multicasting scenario with B = 3 BSs, M = 2 antennas per BS, and \K.i \ = I/C 2 I = l^sl = 10, 
we show the average group-rates of multi-group multicasting using the proposed, the subset search, the broadcast scheme. The 
results of both full BS cooperation and interference coordination scenarios are shown. 


user in the target user subset, users that are elose to the BS (but outside of the target subset) 
may also have the opportunity to suceessfully deeode. Henee, the number of users served is 
often greater than the number of users selected. However, we see from Fig. 8 that there is not a 
significant difference between the two because of the directionality of the signal (i.e., the choice 
of the transmit covariance matrix). Interestingly, this is different from the single-antenna scenario 
where the two numbers may have a significant difference since, without spatial directionality, 
users close to the BS will have a high probability of being served when the target user subset 
includes a cell-edge user. 

B. Multi-Group Multicasting Scenario 

In this section, we consider the multi-group multicasting scenario with B = 3 BSs serving 
collaboratively G = 3 multicast groups, each with 10 users (i.e., |/Ci| = I/C2I = I/C3I = 10). In 
Fig. 9, we show the average group-rate of the relaxed GRB via SCA scheme proposed in Algo¬ 
rithm 4 (c.f. Section IV), and compare it with the subset search algorithm (averaged over only 
300 channel realizations) and the broadcast scheme, where all users are served simultaneously in 
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Fig. 10. For the multi-group multicasting scenario with B = S BSs, M = 2 antennas per BS, and |/Ci| = \K .2 \ = lA^s] = 10, 
we show the average number of selected users per group in the proposed and the subset search algorithms. The results of both 
full BS cooperation and interference coordination scenarios are shown. 


each slot. Both cases with full BS cooperation and interference coordination are considered. In 
the latter case, we assume that each BS serves only one group (namely, BS b serves group g, for 
b = g e {1, 2, 3}) and, thus, {Qg}?,,?,' = Omxm, for all b ^ b' and for all b ^ g. Different from 
the single-group scenario, we can see that user selection in the multi-group scenario is more 
advantageous in the high SNR regime and the gain is much more significant than that in the 
single-group case. This is because, in the multi-group scenario, the performance is interference 
limited at high SNR and, thus, user selection not only can help avoid rate limitations by the 
user with the worst channel but can also help reduce interference between signals intended for 
different groups. More interestingly, the proposed scheme also outperforms the subset search 
algorithm in the multigroup scenario since the latter scheme is more likely to converge towards 
a locally optimal solution in this scenario. These advantages can be observed in both full BS 
cooperation and interference coordination scenarios. Interestingly, the group-rate improvement 
is more significant for the case with only interference coordination. This is because, when BSs 
are not able to fully cooperate, the spatial degrees of freedom are not sufficient to effectively 
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Fig. 11. For the multi-group multicasting scenario with B = 3 BSs, M = 2 antennas per BS, and |A3i| = I/C 2 I = lAi^sl = 10 
non-i.i.d. users, we show the average group-rates of multi-group multicasting using the proposed fair OUS scheme and the 
broadcast scheme. The results of both full BS cooperation and interference coordination scenarios are shown. The minimum 
and 90% rates of the proposed fair OUS scheme is also shown for comparison. 


reduce interference solely through the design of the transmit covariance matrix. Therefore, the 
benefit of reducing interference through user selection is more pronounced in this case. 

In Fig. 10, we show the average number of selected users per group when using the proposed 
algorithm under full BS cooperation and interference coordination. Interestingly, we can see that, 
different from the single-group scenario, the number of selected users is less at high SNR instead 
of at low SNR. This is due to the fact that, at high SNR, the performance is interference limited 
and, thus, the system would benefit more from eliminating users and reducing interference. This 
effect is more evident in the case of full BS cooperation where more spatial degrees of freedom 
are available for signal enhancement and interference avoidance. 

In Fig. 11, we consider the multi-group multicasting scenario with users whose channel vectors 
are non-identically distributed. Again, we have B = 3 BSs serving collaboratively G = 3 
multicast groups, each with 10 users (i.e., |/Ci| = I/C 2 I = I/C 3 I = 10). The users are deployed in 
the same way as that in Figs. 7 and 8 . Again, 6 sets of random user locations are considered, 
each averaged over 100 channel realizations. Each BS serves a group consisting of the closest 10 
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Fig. 12. For the multi-group multicasting scenario with B = 3 BSs, M = 2 antennas per BS, and |A3i| = I/C 2 I = lA^isI = 10 
non-i.i.d. users, we show the average number of selected and served users per group obtained with the proposed fair OUS 
scheme and the broadcast scheme. The results of both full BS cooperation and interference coordination scenarios are shown. 


users. In the figure, we show the average group-rates of the proposed OUS seheme (e.f. Seetion 
IV-B) and the broadeast seheme. The minimum rate among all users and the 90% rate are also 
shown for eomparison. Both cases with full BS cooperation and interference coordination are 
considered. Similar to the i.i.d. case, we can see that the advantages of OUS increase with SNR 
and the gains are much more significant than the single-group scenario. Moreover, a significant 
advantage can still be observed in the case of interference coordination. This is because the 
normalized channel vectors used in the proposed OUS scheme preserves the direction of the 
channel vectors and, thus, is still able to successfully perform interference coordination among 
different BSs. Similar to the single-group scenario, the minimum rate may be smaller than the 
rate achieved in the broadcast scheme if the transmission does not occur over a large number of 
time slots, which is the case in the full BS cooperation scenario. However, the majority of users 
(in fact, over 90% of users) achieve rates that are higher than that of the broadcast scheme. In 
Fig. 12, we show the average number of selected and served users obtained using the proposed 
fair OUS scheme under both full BS cooperation and interference coordination. We can see that 
the number of served users is again more than that of selected target users. However, this effect 
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Fig. 13. Multi-group multicasting scenario with B — 3 BSs, M = 2 antennas per BS, and \K.i \ = I/C 2 I = lA^sl = 20, among 
which only 5 per group feedback their channel vectors. We plot the average group-rate of the cases with and without virtual 
users. 


is less pronounced under full BS cooperation since the signal in this case contains directionality. 

Finally, let us consider the GRB-PF problem as discussed in Section V. Here, we assume that 
there are 20 users per group (i.e., |/Ci| = I/C2I = I/C3I = 20), but only 5 per group feedback 
their channel vectors in each block. Notice that, in addition to the proposed scheme, it is also 
possible to compute the input covariance matrix, the rate, and the user selection assuming that 
only the users who feedback their CSI exist in the network. The latter is referred to as the case 
without virtual users. In the experiments, the proposed scheme is implemented with J = 100 
virtual users. In Fig. 13, we show the average group-rate achieved when the system parameters 
are derived using the above two approaches. We can see that, by considering virtual users, the 
group-rate can be significantly improved, especially at high SNR where user selection is critical. 
The performance of the broadcast scheme where the system parameters are designed by serving 
all feedback users in each block is also plotted for comparison. We can see that, at high SNR, 
where the performance is interference limited, user selection (even without consideration of 
virtual users) can provide significant group-rate improvement. 
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VII. Conclusion 

In this work, the OUS scheme was examined for both single-group and multi-group multicas¬ 
ting scenarios in the physical layer of a multicell multi-antenna wireless system. User selection 
along with application layer erasure coding was proposed to overcome rate limitations caused by 
the worst user in the multicast group. The proposed user selection policies were derived based 
on the optimization and relaxation of a set of user selection variables. For the single group 
scenario, we formulated the problem as a group-rate maximization problem, and proposed an 
efficient user selection policy by performing a convex relaxation and by employing a sequential 
deflation algorithm. For the multi-group scenario, we formulated the problem as a group-rate 
balancing problem and proposed an efficient user selection policy by performing SCA along 
with the sequential deflation algorithm. When only part of the users feedback their instantaneous 
CSI, we further introduced the concept of virtual users to take into consideration the probability 
that non-feedback users are served in each block. The effectiveness of the proposed schemes 
was shown via computer simulations. Interestingly, we observed that user selection is more 
advantageous in the low SNR regime for the single-group scenario, but is more beneficial in the 
high SNR regime for the multi-group scenario. 


We first show that, for S < 


Appendix A 
Proof of Lemma 1 

maxfc log2(l -f Y.b=i Pb\\iib,k\\‘^) 


1 -1 


, (Q, R, {sk}k&!c) is a feasible 


point of (13) if and only if {Q,R,A), where A = {k E 1C : Sk = 1}, is a feasible point of 
(10). Specifically, let (Q, R, {sk}k&K.) be a feasible point of (13) and let A = {k E JC : Sk = 1}. 
We can see that, for any k E A, the constraint in (10b) is equivalent to the constraint in (13b) 
when Sk = 1. Therefore, if (Q, i?, is a feasible point of (13), then (Q,i?, ^) must be 

a feasible point of (10b) as well. On the other hand, let (Q,i?, ^) be a feasible point of (10b) 
and let Sk = I If k E A, and = 0, if /c ^ A. Similarly, the constraints in (13b) are the 
same as those in (10b), for A; G ^ (i.e., for k such that Sk = 1). However, for k' ^ A (i.e., 
for k' such that Sk' = 0), the constraints in (13b) are redundant since log 2 [l -f tr(Qhfc/h|()] -f 
(5“^(1 — Sk') > log 2 [l -I- tr(Qhfch|^)] > R, for all k E A. Therefore, {Q, R, {skjkeK) is a 
feasible point of (13) if (Q,i?, ^) is a feasible point of (10b). Moreover, one can also see that 
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(Q,-R, and (Q,-R, {sfej^gx:) achieve the same objective values in their respective problems 
sinee j^RY^keic^k = Efceyc The lemma follows. 

Appendix B 

Proof of Proposition 1 

We basieally show that Algorithm 4 is a special case of the sueeessive upper-bound mini¬ 
mization (SUM) method in [39]. 

First, notice that problem (21) is the epigraph form of the following problem: 


max mm < - 

9 


mm[rk{{Qe}f=i) + S ^(l-s^)] 


subjeet to 0<Sk <1, Vfc G /C, (20o), and (20d), 
variables: {Qg}^=i,{sfc} 


IkeK- 


(36a) 

(36b) 

(36e) 


Sinee Rg and a are auxiliary variables, we can focus on showing that any limit point of 
{{Qg}g=i 5 is a stationary point of problem (36). Similar to problem (21), problem 

(25) is the epigraph form of 


E 


max mm 

9 


fce/Cg 


^g\^9 


mm 


^fc({Q^}^=il{Q4Ei)+ <i ^(i-Sfc) 


subjeet to 0<Sfc <1, Vfc G /C, (20o), and (20d), 
variables: {QgjJLi, {sfclfcGX- 


(37a) 

(37b) 

(37c) 


The approximation rk{{Q,e}l=i) ~ I {Q4£=i) satisfies rk{{Qi}%i) = rk{{Qi}% 

?^fc({Q£}?=i) > rk{{CR}%i I and 


'9rfc({Qf}F=i 


<9Qs 


0r,({Q,}ti I {Qa.ti 


Qi=Q«,V£ 


RQg 


Q£ = Q£!V£ 


for all g. Therefore, Algorithm 4 is essentially the SUM method [39]. Aeeording to [39, Theorem 
1], any limit point of {{Qg}^=i, {SfcjfcLi} generated by Algorithm 4 is a stationary point of (36). 
Proposition 1 is thus proved. 
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